Nuclear Inelastic X-Ray Scattering of FeO to 48 GPa 
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The partial density of vibrational states has been measured for Fe in compressed FeO (wustite) 
using nuclear resonant inelastic x-ray scattering. Substantial changes have been observed in the 
overall shape of the density of states close to the magnetic transiton around 20 GPa from the 
paramagnetic (low pressure) to the antiferromagnetic (high pressure) state. Our data indicate a 
substantial softening of the aggregate sound velocities far below the transition, starting between 
5 and 10 GPa. This is consistent with recent radial x-ray diffraction measurements of the elastic 
constants in FeO. The results indicate that strong magnetoelastic coupling in FeO is the driving 
force behind the changes in the phonon spectrum of FeO. 
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The study of the electronic and magnetic properties 
of simple transition-metal compounds is an important 
topic in diverse fields ranging from solid-state physics to 
Earth sciences. Iron oxide FeO (wustite) belongs to the 
group of highly correlated transition metal oxides, be- 
ing an archetypal insulating antiferromagnetic material 
at zero temperature. NiO, CoO, MnO fall in the same 
group of materials, which are still not well understood by 
theory. They also present a challenge to experimental- 
ist because the predicted high pressure transformations 
in these materials (metallization, high-spin to low-spin 
transitions) may occur at extreme pressures ( e. g. > 
100 GPa). FeO stands out from this group of materials 
in the sense that it is a possible major constituent of the 
Earth's lower mantle and core, thus its pressure and tem- 
perature dependent properties are very important for the 
understanding of the Earth's interior. Here we present 
a study of the vibrational density of states of isotope- 
enriched Feo.g470 using high-resolution nuclear resonant 
inelastic x-ray scattering. We observe changes in the den- 
sity of states that are consistent with the softening of 
the aggregate shear sound wave velocity of wustite under 
pressure [Q, which we associate with the phase transi- 
tion from the cubic paramagnetic phase to the rhombohe- 
drally distorted antiferromagnetic phase stable at higher 
pressures. We show that the phonon density of states is 
affected by magnetoelastic coupling between phonon and 
spin subsystems in FeO. 

Experiments were performed at the synchrotron beam- 
line SRI-CAT 3ID of the Advanced Photon Source 
(APS), Argonne National Laboratory. The details of the 



diamond cell and the experimental setup are reported 
elsewhere [Q. The sample was ~ 25 /im in diameter 
(loaded without pressure medium), and diamonds with 
flat 400 /xm culets were used. The maximum counting 
rate in the phonon wing ranged from 10 cps at 0.9 GPa 
to 2 cps at 48 GPa. The measured spectra, i. e., count 
rate as a function of the energy difference between the in- 
cident photon energy and the nuclear transition energy, 
I(AE), were converted to phonon DOS profile accord- 
ing to the data analysis procedure described by Hu et al. 
||. Typically, 10-20 DOS spectra (one hour/spectrum) 
at the same pressure were summed together. The pres- 
sure was measured using the ruby fluorescence line and 
a nonhydrostatic pressure scale || . 

Experimental data are presented in Fig.|l|. In Fig.^j we 
compare the phonon DOS of FeO as derived from neu- 
tron scattering || and the partial phonon DOS of Fe at 
0.9 GPa from our measurements. High frequency oxy- 
gen modes are absent from our data, indicating that Fe 
vibrational amplitude is negligible at oxygen vibrational 
frequencies. One apparent feature of the experimental 
data is an inelastic peak at small energy transfers which 
develops at 10 GPa, and persists through almost the en- 
tire pressure range of the present experiment to 48 GPa. 
It is most pronounced near 16 and 20 GPa, which is ex- 
actly the pressure range of the reported transition from 
the paramagnetic rock-salt-type structure to the antifer- 
romagnetic rhombohedral structure . The transition 
is very sensitive to (nonhydrostatic) stress conditions ; 
for example, the splitting of the (111) diffraction line be- 
ing smeared between 10 to 18 GPa if no pressure medium 
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was used, and occuring around 16 GPa § to 17 GPa |] 
under quasihydrostatic conditions. 

The calculated partial density of states weighted by 
the square of the energy E 2 is plotted versus energy in 
Fig.||. In FigJ|, the density of states in the Debye ap- 
proximation would be a straight horizontal line, related 
to the aggregate sound velocity of the material . The 
apparent feature of our data is the pronounced softening 
of the low energy vibrational spectrum. The resolution 
of the present measurements (2.4 meV) does not allow us 
to follow the softening to the low frequency region, where 
most static and ultrasonic measurements are performed. 
However, this behavior suggests large effects on the static 
elastic constants at the Neel transition. 

The temperature dependence of the elastic constants 
of Feo.920 at ambient pressure was investigated in detail 
by Sumino et al. || . Substantial softening of the shear 
constant C44 was found close to the Neel transition. Sim- 
ilar anomalies have also been observed in MnO (softening 
of C44) Jll|. The magnetic structure studies p5|-p7|| by 
neutron scattering have shown that the Mn spins align 
ferromagnetically within a given (111) plane, and these 
planes are stacked antiferromagnetically in the <111> di- 
rection. Similarly, FeO has the same magnetic structure, 
with magnetic moments pointing in the <111> direction. 
At temperatures below the Neel temperature MnO un- 
dergoes a lattice contraction along the antiferromagnetic 
stacking direction <111> pl|. The distortion in FeO 
has different sign, the elongation along the body diago- 
nal. The magnetic system breaks into multiple domains 
(called T domains) corresponding to the antiferromag- 
netic stacking along the four possible <111> directions 
within the crystal. Experimental p8[-p0| and theoretical 
pl| p3[ studies in MnO invoke exchange-striction effects 
to explain the character of the transition. The results of 
these studies can be summarized as follows. The energy 
balance is determined by minimizing the sum of mag- 
netic exchange energy and elastic energy, resulting in a 
rhombohcdral distortion 
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Here A denotes the distortion of trigonaly deformed cube 
corner angles ^7r ± A, C44 is the shear constant, J\ is 
the exchange integral for the nearest neighbors, z\ is the 
number of the nearest neighbors, t\ — —rSlnJi/Sr, and 



S —< Si ■ Sj >f m 



< Si ■ Sj >°„ is a correlation func- 



tion for nearest neighbors with parallel spins (residing 
within the same <111> plane) and antiparallel spins in 
neighboring <111> planes. The numerical value of A is 
about 0.5 degree at 4 K in MnO. Actually, according to 
Bartel, J22j A is proportional to the sublattice magneti- 
zation and is a sensitive measure of the order parameter 
of the magnetic phase. 

However, FeO and CoO have much larger volume 
anomalies below Tjv and their magnetic properties ac- 
cording to Kanamori |24 23] cannot be described by the 



theory developed by Lines and Bartel |22]. In a crys- 
tal field of cubic symmetry, the orbital degeneracies of 
Fe 2+ and Co 2+ are not completely removed, and the 
residual orbital angular momenta contribute to the en- 
ergy balance through spin-orbit coupling and the direct 
effect of orbital magnetic moments, resulting in observ- 
able magnetostriction effects. Substantial contribution 
to the magnetostriction effects in FeO is due to the mag- 
netic anisotropy energy pSfl . Kanamori |2^| calculated 
the equlibrium strain components for FeO resulting from 
magnetostriction. He used the general formulation de- 
rived by Kittel for cubic systems |2(| (all dj 's are cubic 
elastic constants, B\ and B 2 are cubic magnetoelastic 
constatnts) : 



e u = B l ((1/3) - a?) /{C n - C 12 ), 
e tj = -B 2 aiOij/Cte 
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where are direction cosines of the magnetization. 
Kanamori's treatment leads to an elongation along [111] 
diagonal in FeO, which agrees with experimental results 
p7| . From the experiments on FeO at 95 K J2^] one de- 
rives a rhombohedral angle of 59° 32', or ea;j,=0.705 %. 
To estimate the rhombohedral distortion resulting from 
the exchange interaction (Eqjl]) in FeO, we need reliable 
information about the exchange integrals and their de- 
pendence on the lattice parameter, which to our knowl- 
edge is not available at the moment. The experimental 
values of distortion Jl],[| seem to be quite close to the 
Kanamori's calculation. Moreover, rhombohedral distor- 
tion is enhanced at high pressures ^] , indicating increas- 
ing B2 in FeO under pressure. 

The cubic constant B 2 is responsible for the magne- 
toelastic coupling between the phonon and the magnon 
branches, and our results on the enhanced phonon den- 
sity of states are directly related to the magnetoelastic 
coupling in FeO. The expression for phonon dispersion 
including the magnetoelastic coupling follows from the 
equations of motion for the magnetic moment f29|,|30| : 



(^-c 2 n )(c 2 -^ 2 ) 
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Here lj s is the frequency of sound wave, and = 
flgM (/3 + ak 2 ) is the spin wave frequency in the ab- 
sence of magnetoelastic coupling; p is the density of the 
material, Mq is the magnitude of the magnetic moment 
of one sublattice, — gAIo({3 + 27). The parameters 
a, 7, and (3 are related to the exchange interactions and 
magnetic anisotropy in the material |29), g is the gyro- 
magnetic ratio, and A: is a wavevector. 

The calculated dispersion relations for FeO at 28 GPa 
with magnetoelastic coupling included [[Hj are shown in 
Fig. 4a (estimated at T=295 K), corresponding sound ve- 
locities are shown in Fig. 4b. Dispersion relations (Fig. 4a) 
change by less than 1 meV, well within the resolution of 
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neutron inelastic scattering experiment. However, as is 
evident from Fig. 4b, the effect on the sound velocities is 
more pronounced, being almost 20-30 % within the en- 
ergy transfer range up to 5 meV. This agrees reasonably 
well with the enhancement of the density of states which 
we observe in our experiment (Fig. 3). Quantitave agree- 
ment may be sought along the lines of more elaborated 
theoretical models, similar to Ref. p3]. 

At present we do not have enough experimental infor- 
mation and theoretical understanding to better constrain 
the magnetoelastic coupling at T/T^ ~ 1. However, 
as follows from our observations, the effect is most pro- 
nounced close to the Neel transition. Our observations 
support the Mossbauer measurements, where the mag- 
netic transition was observed starting from 8 GPa un- 
der nonhydroctatic conditions. The reason is that uniax- 
ial strain along the body diagonal may induce magnetic 
moments well below the transition point determined un- 
der hydrostatic conditions (17 GPa). This is consistent 
with the notion that this second order phase transition is 
smeared by the external field (uniaxial strain), which is 
proportional to the order parameter of the broken sym- 
metry phase [ f35| . 

In summary, we have observed substantial softening of 
the density of states at energy transfers below 10 meV in 
FeO at pressures close to 15-20 GPa, which persists in its 
antiferromagnetic phase up to 40-48 GPa. We relate the 
observed softening to the effect of the magnetoelastic cou- 
pling in this material. The theoretical estimates [ p5|j22| 
show that both rhombohedral distortion and magnetoe- 
lastic coupling include substantial contributions from the 
magnetic anisotropy energy in FeO, which arises from 
the spin-orbit coupling in the orbitally degenerate ground 
state |24|,^,^7j . Further experimental work is required to 
clarify the effect of exchange-driven magnetoelastic con- 
tribution. 
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FIG. 1. Inelastic part of the signal in FeO as a function of 
pressure (elastic peak is subtracted). The region close to zero 
energy transfer is shown using bold lines. Note enhanced den- 
sity of states close to T/Tjv ~ 1. Ratio T/Tjv was estimated 
using Tjv=198 K at P=0.1 MPa, and data from Ref.[7]. 
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FIG. 2. Dispersion of phonons and magnons in FeO from 
neutron inelastic scattering experiments. Upper part shows 
comparison of phonon DOS calculated from neutron measure- 
ments (solid line) with our partial DOS for iron at 0.9 GPa 
(grey circles), 300 K. Also shown are measured and calculated 
phonon (circles and solid lines) and magnon (diamonds and 
dashed lines) branches from Ref. [5]. 
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FIG. 3. g(E)/E 2 derived from partial DOS in FeO as a 
function of pressure. For the Debye model, the lower energy 
part should be a horizontal straight line. 
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FIG. 4. Model for the magnetoelastic coupling in FeO. 
a) The interacting transverse phonon and magnon branches. 
Nonintcracting bare frequencies are shown with dashed lines, 
the dispersion branch for magnons is calculated accord- 
ing to [31], the phonon dispersion was approximated by 
E(Q) = 2Qov s /nsin((ir/2)(Q/Qo)) using the sound velocity 
3.4 km/sec at Q=0. b) Calculated sound velocity including 
magnetoelastic coupling (solid line) and without magnetoe- 
lastic coupling (dash-dotted line). 
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